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Abstract. A long range Weak Equivalence Principle (WEP) violating force between Dark Matter (DM) particles, mediated 
by an ultralight scalar, is tightly constrained by galactic dynamics and large scale structure formation. We examine the 
implications of such a "dark force" for several terrestrial experiments, including Eotvos tests of the WEP and direct-detection 
DM searches. The presence of a dark force implies a non-vanishing effect in Eotvos tests that could be probed by current 
and future experiments depending on the DM model. For scalar singlet DM scenarios, a dark force of astrophysically 
relevant magnitude is ruled out in large regions of parameter space by the DM relic density and WEP constraints. WEP 
tests also imply constraints on the Higgs-exchange contributions to the spin-independent (SI) DM-nucleus direct detection 
cross-section. For WIMP scenarios, these considerations constrain Higgs-exchange contributions to the SI cross-section to be 
subleading compared to gauge-boson mediated contributions. The combination of observations from galactic dynamics, large 
scale structure formation, Eotvos experiments, DM-direct-detection experiments, and colliders can further constrain the size 
of new long range forces in the dark sector. 
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INTRODUCTION 

Despite the large body of evidence for Dark Matter (DM) from the galactic rotation curves [HE), acoustic oscillations 
in the cosmic microwave background J3]|4l|5]|6l, large scale structure formation G][8), and gravitational lensing [|9l [T0ll 
almost nothing is known about the properties of DM and a great deal of experimental effort is under way to unravel its 
nature. Ground based direct detection experiments lfTTl[T2l put limits on the DM mass and the strength of its interaction 
with baryonic matter from observations of recoiling nuclei. Experiments lfT3l [141 [T5l [T6l studying cosmic rays from 
the galactic halo have recently seen indications of an electron/positron excess, which could be interpreted as evidence 
for DM annihilation, and can constrain the DM mass and interactions. 

Another set of experiments are devoted to question of whether DM violates the Weak Equivalence Principle (WEP). 
Such an apparent violation of the WEP can occur in the presence of a new long range scalar force coupled to DM 
and is the focus of this work. In particular, we consider the possibility of an ultralight scalar with mass m§ < 10~ 25 
eV coupling directly to DM with intergalactic range. There are a few hints for the presence of a dark force from the 
higher than predicted supercluster densities IfTTl . and voids llT8l 1191 , and galaxy clusters [20] ( for a summary see 
ETI I22T ). Strong constraints on a dark force are derived from observations of DM dynamics in the tidal stream of 
the Sagittarius dwarf galaxy ETI |231 . which indicate a force with strength less that 20% of gravity. However, new 
observational systematic errors have been recently discovered 1 24 1 that could require a revision of this result, perhaps 
allowing for a stronger dark force. A more recent analysis (23 considers the effect of a dark force on the evolution of 
density perturbations and the resulting impact on the CMB spectrum with similar constraints. 

Many models that contain the interaction of an ultralight scalar with DM lf26l [27l [28l |29l [30l [3T1 fTSl [32l [33l [34l 
[19] [35 1 have been proposed to explain features in the DM distribution and explore the possibility of DM-quintessence 
interactions. More recently, work with non-universal scalar-tensor theories of gravity with the Abnormally Weighting 
Energy (AWE) Hypothesis ll36l l37l also invoke couplings of an ultralight scalar to the dark sector as a way of 
explaining the observed cosmic acceleration even in the absence of a dark energy fluid. Constraints on such scenarios 
from big bang nucleosynthesis have also been studied [38 1. 

If the DM has interactions with the Standard Model (SM), the presence of a scalar dark force will be communicated 
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FIGURE 1. Two-loop diagrams in WIMP DM models that generate the operators ffj in Eq. |7J( which generate couplings of the 
ultralight scalar to the SM fermions after electroweak symmetry breaking. 



to the SM via virtual DM. This effect can be exploited to connect signatures of dark forces in astrophysics and cos- 
mology with terrestrial experiments such as the laboratory Eotvos tests of the WEP, DM-direct-detection experiments, 
and in some cases even at colliders. This connection has been studied in 1221 and 1391 and more recently an extensive 
analysis was done in BUI and is the focus of this note. 



WEP VIOLATION PHENOMENOLOGY 



We first discuss some basic phenomenology of new long range Yukawa forces coupling to ordinary matter or DM. The 
potential V between two bodies of mass M, and M s , in the presence of a new long range Yukawa force mediated by a 
scalar of mass m^, can be parameterized as 
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where the first term in the potential is the usual Newtonian gravitational potential and the the parameter Gfe is a 
product of the charge to mass ratios 2i,. v /M,, s of the two bodies under the new Yukawa force. The factors of fy zS appear 
for dimensional reasons [40 1 and are 1 for fermionic test objects and l/(2Mj. s ) for scalar test objects. The Eotvos 
parameter 77^ which measures the relative acceleration between two test objects in the presence of a gravitational 
source s is given by 



m -ai\ 



1 



4nG 



G2I2 



Mi 



M 2 



Q.A 



(2) 



where we ignored the small mass < 10 25 eV for the ranges r <C being probed. The current bounds ATI on 

the Eotvos parameters TJ^ e ' Tl and T?^' Tl which measure the relative acceleration of macroscopic samples of Beryllium 
and Titanium toward the center of the Earth and galactic DM respectively are 



^Beji < (0.3 ±1.8) X 10" 
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The Apollo (LLR)||42] and Microscope 1 43 1 collaborations are expected to be sensitive to values of Tj^ e ' Tl as small as 
10~ 14 and 10~ 15 respectively. The MiniSTEP[44| proposal, if approved, can be sensitive to values of rj^ e ' Tl as small 
as 10~ 18 . 

We parameterize the coupling of the ultralight scalar to DM by 
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FIGURE 2. An estimate of the allowed region in the (j] OM£ ,/3) parameter space for minimal WIMP DM models. The curves 
in the figures give an estimate of for rj DME for a given value of j8 from the two loop diagrams in Fig. [T] The shaded region is 
unlikely for typical WIMP models. Using the observational constraint j6 < 0.2, the allowed region is further restricted to the left 
of the vertical line. The estimates in the above figures for IJ™'? , for /3 < 0.2, are far below the current experimental bounds 

^dm T ' < 10~ 5 , We < 10~ 13 . An improvement of about five orders of magnitude would be required in Eotvos experiments to fully 
probe the allowed parameter space for jS = 0.2 by measuring 7]^ e ' Tl . This is within reach of the MiniSTEP proposal|44 | . 



where % denotes the DM field and g x is the coupling of DM to the ultralight scalar. For simplicity we assume that 
fermionic DM lives in vector-like representations of the electroweak gauge group so that the interaction is gauge 
invariant. For chiral DM the interaction with the ultralight scalar can arise only from higher dimension operators. 
From observations in astrophysics 12T1I231 and cosmology [25] , the most recent bound on the charge to mass ratio of 
DM under the new dark force parameterized by /3 is given as 

V47T M x 

where Mp denotes the Planck Mass. This bound corresponds to a dark force restricted to be less than about 20% of 
gravity. 



SCALAR DARK FORCES AND TERRESTRIAL EXPERIMENTS 

If the DM has interactions with the SM then the bounds on the parameters fj^^ 1 and j3, obtained from laboratory 
WEP tests and astrophysical and cosmological observations respectively, will in general be correlated. As a simple 
example consider minimal WIMP DM models l45l of the type 

_g = \x(U? + Mo)X, fermionic DM, 

\c{D^D^x~cMlxh~V{x,H), scalarDM, 

where a single DM particle % is added to the SM with non-zero electroweak interactions. Relic density considerations 
typically require the DM mass to be in the TeV range in such models. We point out that for non-zero hypercharge, 
such minimal WIMP DM models are typically ruled out RBI by direct detection experiments. Here we consider such 
minimal DM models, even with non-zero hypercharge, only as illustrative examples keeping in mind that such DM 
could be part of a non-minimal extension which avoids the direct detection bounds. If the DM couples directly to a 
dark force as in Eq.Q, then this dark force will couple to the SM quarks and leptons via two diagrams as shown in 
Fig.[T] In particular, before electroweak symmetry breaking these loop diagrams generate the dimension five operators 
of the type 

6* = <$>Q L £rfC^u R , 0$ = <j>Q L HC%d R , 0? = <j>LLHC?e R , (7) 
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FIGURE 3. DM-induced coupling of (j) to SM fermions at one loop in the presence of additional squark and slepton like degrees 
of freedom. 

where the coefficients C~j are finite and computable from the loop diagrams. After electroweak symmetry breaking, 
the Higgs field gets a vacuum expectation value so that these operators generate a coupling between the ultralight 
scalar and the SM fermions. Note that the coupling of the ultralight scalar to the SM fermions must come from 
induced higher dimension operators since no gauge invariant renormalizable couplings exist between SM fermions 
and a gauge singlet scalar before electroweak symmetry breaking. From the induced coupling of the ultralight scalar 
to the SM model fermions one can obtain order of magnitude estimates [40] for the Eotvos parameters in terms of the 
parameter /3. For minimal WIMP models of the type in Eq.fjfrji, the lower bound on the Eotvos parameters fj^ 5 ^ 1 as a 
function of j3 is plotted in Fig. [2] For j3 = 0.2, marked by the vertical lines in Fig. [2] we see that the lower bounds for 
typical WIMP DM models are n* e ' Tl > 4 x 10~ 12 and n Be - Tl > 4 x 10~ 18 . These lower bounds are far below the current 
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experimental upper bounds shown in Eq.lpk. An improvement of about five to seven orders of magnitude in Eotvos 
experiments would be required in order to probe these lower bounds for WIMP DM models. The MiniSTEP [44] 
experiment, which is currently under study, is expected to reach a sensitivity for 77^ e ' Tl of about 10~ 18 and might be 
able to probe the lower bounds. If an effect is detected in J^'^ 1 far above the lower bounds in Fig.Q it would suggest 
the possibility that the coupling of to the SM fermions might have a different origin. Another possibility that might 
explain an effect above the expectation in Fig. [2] would be a stronger induced coupling of (j) to ordinary matter arising 
from non-minimal DM models; for example a one loop coupling of (j) to ordinary matter (see Fig. [3]) in the presence of 
additional squark degrees of freedom |40l . The main point is that given a DM model one can estimate the expected size 
for the Eotvos parameters r\ E DM for a given value of /3 thus correlating observations in astrophysics and cosmology 
with laboratory WEP tests. 

Another mechanism by which the ultralight scalar can couple to SM fermions is by mixing with the Higgs. While the 
ultralight scalar cannot have any renormalizable couplings to the SM fermions before electroweak symmetry breaking, 
it can have direct couplings to the SM Higgs. After electroweak symmetry breaking, the ultralight scalar will mix with 
the SM Higgs and in turn couple to the SM fermions. The size of this mixing angle will be tightly constrained by 
Eotvos tests since such an induced coupling to the SM fermions will lead to a violation of the WEP. In what follows 
we will refer to the ultralight scalar gauge eigenstate as S and the mass eigenstate after electroweak symmetry breaking 
as <p . The renormalizable couplings of the ultralight scalar S to the SM Higgs doublet are encoded in the potential 



We assume for simplicity that S does not acquire a vaccum expectation value. After electroweak symmetry breaking 
the H T HS interaction induces mixing between the Higgs boson and the ultralight scalar. In unitary gauge the neutral 
component of the Higgs doublet H is given by 



V(H,S) 




(8) 



•o 



v + h 




(9) 



/ \ / \ ' \ H T V''H 

/ * ' \ ' \ ' \ 

' 1 ' 1 ' « ' I 



I » i "•• + ; x i l ± ; x i l | ; x | L. i 

1 ' , I I I I ' 

\ ' \ / V / \ ' 

h;..h h/\h 



FIGURE 4. One loop diagrams which which contribute to terms in the effective potential V(H,S) that are linear in the 5 field. 
The first two diagrams are UV divergent and contribute to the renormalization of the 5-tadpole and the coupling Si respectively. 
The remaining diagrams mix into higher dimensional operators and give a finite contribution to Higgs-ultralight-scalar mixing after 
electroweak symmetry breaking. 



and the mass terms in the potential become 



V mass = \{nlh 2 + n 2 s S 2 + nl s hS), (10) 



where 
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Mft = ^-, Vs = *2 + -^-, MfcS = SlV. (11) 
The mass eigenstates h± in terms of S and h can be written in terms of a mixing angle 9 as 
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with corresponding masses 
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We identify h + and m + with the physical Higgs boson field and mass respectively and the physical ultralight scalar 
field and mass are identified as 

= /z_, m ( j l = m_. (14) 

Now if the DM couples to the ultralight scalar as in Eq.|4]i, the mixing angle 9 could receive additional contributions 
from virtual DM effects. As a concrete example consider the DM model where the DM % is either a real singlet or real 
triplet scalar under the electroweak gauge group with potential 

V(H,S,X) = V(//,5)+iM V + ^Z 4 + «2//W+^ 2 5 + A^ 2 5 2 , (15) 

where we have imposed the Zj_ symmetry % — ► — % for the DM field to ensure its stability. The experimental constraints 
on this DM model in the absence of a dark force were studied in |46 47j. The parameter ai determines the spin- 
independent cross-section for DM direct detection 

2 2 2 2 

a * N - „rJ 1 V2 4 - M x = M o + "2v 2 , (16) 




FIGURE 5. Upper bound on 02 in the singlet (upper curve) and real triplet (lower curve) scalar DM models as a function of 
T7^ e ' Tl /j3 2 . The vertical black line in the left panel corresponds to the current bound of r]^ e ' T ' < (0.3 ±1.8) x 10~ 13 . This vertical 
black line will move to the left with further improvements in Eotvos experiments as indicated by the arrow pointing to the left. The 
right panel shows the region closer to the expected future bounds from the MiniSTEP experiment and the vertical black line here 
corresponds to the expected sensitivity of MINISTEP of rj^ eJi = 1(T 18 . We have used m h = 120 GeV and /3 = 0.2. 



where M% is the physical DM mass after electroweak symmetry breaking. For the singlet DM model, 02 also determines 
the relic density as a function of the Higgs mass. We now show that the presence of a dark force can imply an upper 
bound on the magnitude of which can translate into an upper bound on the spin-independent direct detection cross- 
section. For the singlet DM model, the implied upper bound on the magnitude of can be used to rule out dark forces 
in large regions of parameter space from the DM relic density constraint. 

It was shown in [40 1 that the mixing angle between the ultralight scalar and the Higgs is of the form 

a 2 v 3 8 v 

sin0 wtan0 k-^- y /3 + ^y, (17) 

M P ml rrr h 

where K = 1 , 3 for the real singlet and real triplet DM respectively. The second term on the RHS above comes from 
the renormalizable term H^HS term after electroweak symmetry breaking. The first term on the RHS in Eq.( 17 1 is 



the result of summing the contributions of higher dimension operators of the form (H'H)"S, for n > 2, induced via 
virtual DM as shown in Figj4] at one loop. The first two diagrams in Fig. |4] are UV divergent and just contribute 
to renormalization of the S-tadpole and the 8\ coupling respectively. The remainin g di agrams are finite and induce 
operators of the form (H^H) n S which give rise to the first term on the RHS of Eq.( 17 1 after electroweak symmetry 
breaking. For more details see l40l . 

As mentioned earlier, Eotvos experiments will constrain the size of sin 9 since it leads to WEP violation for ordinary 
matter and will generate a non-zero rjf e 'J^. Furthermore, note that the RHS of Eq.( 17 1 involves the parameter /3 which 



is constrained from cosmology and astrophysics as in Eq.Q. Except in slices oTparameter space where there are 



intricate cancellations between the two terms on the RHS of Eq.( 17 1, which is unlikely given that 81 and 02 are not 
related in any way, one can use the Eotvos tests to put an upper bound |40| on the magnitude of «2 as shown in Fig. 5 
Figure |5] shows the upper bound on the magnitude of 02 as a function of the ratio e ' Tl //3 2 . In both panels of Fig. 5 
the upper curve is for the singlet scalar DM model and the lower curve is for the real triplet scalar DM model. The 
vertical black line on the right end of the left panel in Fig. pi corresponds to the values rjf^ r ' = (0.3 ± 1.8) x 10~ 13 
and j3 = 0.2 which are the current bounds. The vertical blacK line in the right panel in Fig fBl corresponds to the value 
Tjg e ' Ti — 10~ 18 , which is the expected sensitivity for the MINISTEP proposal, and j5 = 0.2. For singlet scalar DM, the 
upper bound on «2 from Fig. [5] can be compared with the required value of a2 l47l for a given Higgs mass in order 
to get the right DM relic density and thus further constrain dark forces. For example, for a Higgs mass of 120 GeV 
and singlet scalar DM mass of 20 GeV, the required value for 02 is about 0.15. In this case a dark force with strength 
corresponding /3 = 0.2 is ruled out by the current bound on TJ^ e Tl shown in Eq.([3j). This can be since in the left panel 
of Fig. [5] which requires 02 < 0.045 as seen from the intersection of the upper curve with the vertical line on the right. 




FIGURE 6. Upper bounds on the spin-independent direct detection cross-section for the scalar singlet DM model as a function 
of r/^ e ' Tl / j8 2 . The vertical line in the left panel corresponds to the current bound of ri^ e Tl < (0.3 ± 1.8) x 10~ 13 . This vertical line 
will move to the left with further improvements in Eotvos experiments as indicated by the arrow pointing to the left. The right panel 
shows the region closer to the expected future bounds from the MiniSTEP experiment and the vertical line here corresponds to the 
expected sensitivity of MINISTEP of rj^ e - Tl = 10~ 18 . We have used j3 = 0.2 and the tree curves in each panel correspond to Higgs 
mass values of m/, = 120, 130, 140 GeV as indicated. 

From Eq.tfTo), the upper bound on the magnitude of 02 also translates into an upper bound on the spin-independent 
cross-section as shown in Fig.|6]for the singlet DM model. For brevity we have not shown the corresponding plots [40] 
for the real triplet DM model. Figure [6] shows how the spin-independent DM direct detection cross-section, Eotvos 
tests, and observations in astrophysics and cosmology can be brought together to constrain dark forces. As seen in 
Fig.pl smaller values of j3 for a given value of J]^ e Tl , lead to a weaker bound allowing for larger values of the spin- 
independent direct detection cross-section. As seen from the left panel in Fig.|6| for a DM mass of around 100 GeV, 
the upper bounds on the spin-independent cross-section are within reach of current or future sensitivities of direct 
detection experiments. For a more detailed discussion see [40] . Depending on the details of the DM model one can 
also perform a similar analysis to correlate observations in WEP tests and astrophysics and cosmology with collider 
signals. One such example for a real scalar triplet as part of multicomponent DM was studied in [40]. 

As long as the DM is not sterile, the presence of a WEP violating scalar force in the dark sector will be communi- 
cated to the ordinary matter sector via virtual DM effects. Depending on the DM model, this effect can be exploited 
to constrain dark forces weaker than gravity through correlated observations in astrophysics, cosmology, laboratory 
WEP tests, DM-direct-detection experiments, and collider signals. 
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